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The use of non-orthogonal, or conic, mounts for steerable antennas in- 
troduces some control system design problems not present in the more 
conventional orthogonal mounts. These problems result from both the 
geometrical and the mechanical cross-coupling which occurs between the 
two non-orthogonal axes of motion. 

This paper presents a general analysis and design of the control system 
for the open cassegrain antenna which can be readily applied to other non- 
orthogonal antenna structures. The form of the feedback controller for ap- 
proximately non-interacting control of each axis is developed. Also de- 
scribed is a supplementary control strategy for providing tracking near the 
zenith region without excessively high slewing rates. 

A computer simulation o1 the system has verified the basic control strategy 
for non-orthogonal mounts and established the feasabilily of operating 
compact antenna structures such as the open cassegrain design under 
severe wind conditions without a radome. 

I. GENERAL SYSTEM DESCRIPTION 

While the general control system design methods developed in this 
study apply to any antenna mount using two non-orthogonal axes, the 
specific structure considered throughout this paper is the slant -mounted 
open-cassegrain antenna. 1 In this structure, the antenna beam tracks 
the target by controlled rotational motion about the inclined and the 
vertical axes shown in Fig. 1. While motion about the vertical axis pro- 
duces true azimuth motion of the beam, the motion about the inclined 
axis generates a combination of azimuth and elevation motion of the 
beam. Further cross-coupling between the azimuth and elevation track- 
ing channels is introduced by unavoidable mechanical coupling of 
motion between the two drive axes. 

Fig. 2 is a general block diagram of the tracking control system. The 
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Fig. 1 — Open cassegrain antenna with two-axis conic mount structure 
showing simplified subreflector structure. 

pointing error is resolved for convenience into the standard azimuth and 
elevation angle errors (this is a conceptual, not physical, portion of the 
system). In the tracking of active repeater satellites, the error signals 
are derived from the waveguide mode detector receiving the satellite 
beacon signal. 2 The horizontal and vertical error signals, e/, and e„ , at 
the output of the error detector are related to the azimuth and elevation 
errors by 



€„ = K h (cos E)(A r - A) 
e v = K v (E r - E) 



(1) 
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where K h and K v are the detector gains, A r and E T are the reference 
azimuth and elevation angles, respectively, of the satellite, and A and 
E are the controlled azimuth and elevation angles of the beam axis 
(electrical boresight) of the antenna. 

All of the equations of motion of the antenna drive system and phys- 
ical structure (considered in detail in Section III) can be represented 
here by the single nonlinear differential equation, 



W = F(w ' u) 



(2) 



where W is the state vector containing as components those variables 
necessary to adequately represent the dynamics of the antenna system, 
and F is the vector-valued function relating the time derivative of the 
state to itself and to the control vector, u (with components U\ and u 2 ). 
From the error signals (1), as well as from feedback signals derived 
from the components of the state vector, W, the controller must gen- 
erate the two antenna drive signals, U\ and u 2 , which control the antenna 
angles, V and /, in such a way as to reduce the tracking error and keep 
the antenna beam automatically "locked-on" to the satellite. Because of 
the complex, nonlinear, multivariate nature of this system, the design of 
the controller cannot be achieved by conventional analytic design meth- 
ods. A preliminary controller design, based on a linear approximation of 
the system dynamics, together with the appropriate coordinate trans- 
formations and supervisory control logic, is discussed in Section IV. The 
final design will be evolved from this preliminary design through an 
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Fig. 2 — Preliminary block diagram for antenna control system. 
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accurate computer simulation of the overall system, its inputs, and its 
environment. 

II. CONIC MOUNT CHARACTERISTICS 

To gain an initial understanding of the tracking requirements on this 
system, the basic transformations between the antenna angles (I,V) 
and the tracking angles (E,A) are needed. Some of this is similar to a 
study on conic mounts by Norton, 3 and his notation is used here. 

Fig. 3 shows the geometry of the conic mount, in particular the struc- 
tural design angles (a,/3), the inclined- and vertical-axis angles (I,V), 
and the elevation and azimuth angles (E,A). 



(inclined axis) 

z 




Fig. 3 — Conic mount geometry. 
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Since V produces only azimuth motion, the only coordinate conver- 
sions needed are those giving the azimuth angle A and the elevation 
angle E produced by /. Fig. 3 indicates that these angles are related by 
rotation about the y-axis through an angle (90° — a), which corresponds 
to the rectangular coordinate transformation: 
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(3) 



For a unit radius, the spherical coordinate equivalents of these rec- 
tangular coordinates are, 



cos E cos A 
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(4) 



Multiplying out the right-hand side of (4), the coordinate conversions 
can be expressed as 



E = sin l (sin a cos — cos a sin /3 cos /) 

. -i T — sin /3 sin / 

A = tan -. i 

_cos a cos p 4- sin a sm /3 cos 

A = V + A. 



7 



(5) 



In order to have complete coverage of the zenith region, it is neces- 
sary that E = 90° when J = 180°. From (5) this occurs if and only if 
a + /3 = 90°. All relationships from here on assume this zenith condi- 
tion. In particular, since a -+- = 90°, let us define, 



a 2 + b 2 = 1 



(6) 



a = sin a — cos /3 
b = sin /3 = cos a 

so the elevation-azimuth expressions (5) reduce to* 

* The antenna mount for the open cassegrain design 1 has a = 42.5° and /3 = 
47.5° which gives a —5° to 90° range in E. Therefore sin a = cos & = 0.67559, and 
sin = cos a = 0.73728, but for simplicity of notation, and somewhat more gen- 
erality, we continue to use a and b. 
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E = sin -1 (a 2 - b s cos I) 



A = tan" 1 [-1 tang)] 



(7) 



A = V + A. 

The inverse of these expressions gives the required antenna mount 
angles, I and V, to produce a given azimuth and elevation of the beam 

axis: 



/ = cos 



( a - sin E\ 
V 1 - « 2 / 

V = A - A (8) 

m , -i/l , / l - 2a 2 + sin E \ 
A = -cr(/)tan ^- f j _ sin S J 



where we define 



a (I) =sgn(sin/) - T^fr - »> 



These relationships between the angles 7, A = A - F, and # are 
plotted in rectangular form in Fig. 4, and in polar form in Fig. 5. It is 
apparent from these figures that there are two pairs of drive angles 
(I,V) corresponding to every pair of tracking angles (E,A).* However, 
in continuous tracking of a moving target, the choice of which pair (I,V) 
to use in pointing the antenna to the given (E,A ) is arbitrary only at 
the beginning of the track, since instantaneous switchover to the oppo- 
site pair is not possible. To switch from one pair to the other may be 
necessary or desirable in certain applications (see Section 4.2), but it 
can be accomplished only by moving the antenna boresight off the 
target for some finite period of time. The exception to this is the unique 
case of the target track which passes precisely through the zenith, at 
which instant the two (I,V) pairs coincide, so that instantaneous 
switchover can be made. 

Finally, the coordinate conversions between the conic mount angular 
velocities (j,V) and the tracking angular velocities (E,A ) are of interest 
in determining the control system requirements. Using (6), (8), and (9) 
we get 

* For example, for the trackiug angles (E = 40°, A = 80°), there are the two 
equivalent pairs of drive angles (7i = 110°, Vi = 145°) and (ft = 250°, Vt - 15 ). 
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Fig. 4 — Resolution of inclined axis angle into elevation and azimuth compo 
nents for antenna parameters: a = 42.5°, fi = 47.5°. 



/ = a(l) 1 - 



On- 



1 + sin E 



>E 



A = - — 5-^~ p; J + tt5 ( cos - 

I + a — o cos / \ 2 

7 = i - A 



(10) 



where 5 (cos 1/2) is the Dirac delta function of (cos 1/2), representing 
the derivative of the step discontinuity in A which occurs whenever 
cos (1/2) = 0, i.e., whenever I = ± (2n + l)7r, n = 0,1,2, ••• (see 
Fig. 4). 

To predict the tracking requirements for the two-axis conic mount, 
it is clear from (8) and (10) that the elevation and azimuth angles 
(E,A), and rates (k,A) are needed. In the tracking of communication 
satellites, these are functions of the orbit parameters of the satellites, 
the locations of the tracking station, and the time reference chosen. 
Computer routines for providing such data are available. 4 For the pur- 
pose of analysis and preliminary design, however, it is desirable to have 
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analytic expressions for the azimuth and elevation angles and rates for 
circular orbits. These expressions are derived in Appendix A. 

III. DYNAMICAL MODEL OF ANTENNA AND DRIVE MOTORS 

As an essential step in the design of the control system, we consider 
now the mathematical model of the antenna mechanical system and the 
drive motors. Although the model must realistically represent the out- 
put angle response (I,V) to input signals, (wi , u 2 ), some simplification 
will be made to ease computer simulation. 

The model of the antenna mechanical system, based on a study by 
Coyne, 5 which was considered to adequately represent the essential 
dynamics, is illustrated in Fig. G. 




J -MASS MOMENTS OF INERTIA 
R -MECHANICAL FRICTION 
L -SPRING CONSTANTS 
T a -APPLIED TORQUES 
$ -SHAFT ANGLES 



BASE 
Fig. (i — Model of antenna mechanical system. 



(13) 
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Using the notation in Fig. 6, the equations of motion are 

T a w = J B A + Li(4 - I) + tfi(0i - i) + RA (ii) 

Li(0i - /) + ^1(^1 - i) = JuJ - m , 
for the reflector structure, and 

7'a (2, - JJ* ~ LAV - B s ) - R,(V - 0,) + RA (12) 

-U(V - 6,) - R 2 (V - d,) = ,h, z -V - aT n (,> + m, , 

for the base structure, where 

m , = — b.J Uz {y cos I — IV sin /) + aJ Uz V, 

m 2 = b 2 Ji vv sin 7(7 sin I + 7 F cos I) 

+ 6 cos 7yix*b(7 cos I - IV sin /) - J lxz (aV - I)]. 

The product mass moments of inertia of the reflector section are defined 
with respect to the x-y-z coordinate system. The J> z >,> is measured with 
respect to the vertical (2 ) axis. 

Due to the large speed variations required in tracking, and because a 
stiff drive is needed to cope with disturbance torques, it is expected that 
hydraulic transmissions similar to the units used to drive the Andover 
horn-reflector antenna 6 will be employed as the drive units. In addition, 
direct gearing is assumed. 

The differential equation which describes the hydraulic drive unit'' is 

K u u = Kj m + K L P + K C P (14) 

where u, d m , and P represent the drive signal from the controller, motor 
shaft angle, and hydraulic pressure. The torque delivered by the hy- 
draulic unit is proportional to the pressure, or 

T D = K T P. (15) 

Equation (14) is valid provided the pressure, P, is less than the maxi- 
mum pressure, P max , which is allowed in the transmission. This condi- 
tion is shown in Fig. 7 where the lines AB and CD represent (14) with 
u = ±iw and AD and BC represent P = ±P, lia x • The line EOF repre- 
sents the operating line or static load line of the hydraulic transmission. 
To illustrate the dynamic operating of the hydraulic transmission, sup- 
pose the operating point is at G with u = h and a large change in drive 
signal to u = h occurs. The motor velocity cannot change instantane- 
ously, so that the pressure in the transmission unit becomes larger. If 
the drive signal change is large enough, the P may exceed P mnx and the 
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Fig. 7 — Dynamic operation of hydraulic transmission. 



new operating point jumps to point H in Fig. 7. At this point, a relief 
valve is actuated which maintains the pressure at P mnx and the torque 
at K T P max until the shaft velocity increases to the value at point /. 
Then, (14) is again valid with u = l\ and the shaft velocity will increase 
to the value at ./.If the drive signal change is such that P < P max , the 
dynamic path is similar to KIM. As a result, in the system simulation, 
the value of P must be monitored and, if it exceeds P max , we must set 
P = Pmnx > i-e., a limiting function. 

A counter-torque arrangement of the motors was assumed to eliminate 
hysteresis effects in the drive systems. 6 This permits the linear gear 
train equations, 



T a = N G T D , 

6 m = N G d , 



(16) 



to apply, where T„ , 6, and N a are the torque applied to each axis in Fig. 
6, the shaft angle at the gearbox output, and the gear ratio, respectively. 
The mathematical model for the computer simulation of the antenna 
structure and the drive motors, represented by (11) through (16) and 
Fig. 7, is shown in detailed block diagram form in Figs. 8 and 9. 
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IV. CONTROLLER DESIGN 

4.1 Approximately Non-Interacting Control 

In order to design the controller for approximately non-interacting 
control of each channel, it is necessary first to represent the drive motors 
and antenna structure in terms of a linear system which is a good ap- 
proximation of the actual system in the small error-signal case and nor- 
mal operating conditions. We assume first that the drive motors are 
unsaturated, i.e. P < P, uax . Next we assume that the torsional spring 
constants of the reflector and base structures are sufficiently large that 
0i (0 - /(0 and 2 (0 = 7(0- 

Although these assumptions linearize and simplify the drive and self- 
coupling portion of each channel in Fig. 8, some major complexities and 
nonlinearities of the system remain in the cross-coupling portion shown 
in Fig. 9. There are, unfortunately, no reasonable or standard assump- 
tions to apply to this portion, only educated guesses. After examining 
the relative magnitude of the various terms in Fig. 9 for typical satellite 
tracks, the simplified linear system representation shown in Fig. 10 was 
chosen for the preliminary design study. 
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Fig. 9 — Cross coupling portion of Fig. 8. 
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Fig. 10 — Simplified linear system representation for deriving basic non- 
interacting controller. 

Fig. 10 also shows the elements of the basic controller design: the error 
signal conversion unit and the feedback compensation units. 

The channel error signals, which are conceptually shown in Fig. 10 as 
the difference between the desired, or reference, angle and the output 
angle for each channel, are obtained as physical signals from the con- 
version of the autotrack error signals, given in (1). This conversion, de- 
rived in Appendix B under the assumption of small errors, is given by, 



62 - [db cs ° 2 fl €v + [ Kkb * sin2 7 ( 1 + fl2 cot2 J" e * - 



(17; 



where K v and Kh are the detector gains in (1 ), and the constants a and b 
are defined in (6). 

The controller outputs, u x and Ma , are derived from these error signals 
plus additional compensation through the feedback networks B(s) and 
C(s). The purpose of these networks is to improve the response to errors 
in the same channel, and to eliminate the undesired response to errors 
in the opposite channel. For the development of this non-interacting 
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design it is convenient to use matrix notation. We define the input, out- 
put and control vectors, 



KvJ' * = \v 



)• -(;> 



The overall input-output relationship is denoted by 

«(•) = D<e r (s) (18) 

where, for non-interacting channels, the transfer matrix must be diag- 
onal. Therefore, we require D to have the form 



"l 



" i D i s) w * ,i: 



From the block diagram of Fig. 10, we obtain the intermediate rela- 
tionships, 

Gie(s) = Fu(s) 



(20) 

U(s) = tQr(s) ~ C<q(s), 



where 



r 'l + Bi(«) CM , 

[ cm i + bmJ' 



( Fl °) 

\aF t Ft/' 



and (21) 

G _ /V + (ffi + h)s -g\s 2 \ 

\ Qtf* + ahtf g 2 s 2 + (03 + h)s)' 

Then, from (18) and (20), it follows that the controller matrix, 0, is 
given by 

C = D' 1 - F~ l G. (22) 

Upon substitution of the matrices in (19) and (21) into (22), we obtain 
the required transfer functions of the controller units: 

*iW = TTTx ~ ' " IF fo* S + (0i + *0«1 
Di(s) /'i 

B 2 (s) = fr-T-r - 1 - w Kg* + a^ 4 )s 2 + (03 + Aj)s] 

* W * (23) 

d(«) - i W) 

C 2 (s) = =- Ua0o - 0s) s 2 + agis] . 
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The only design objective incorporated in these controller functions is 
that of non-interacting channel control. This has fixed the cross-coupling 
controller units Ci(s) and C 2 (s) in terms of the linear plant parameters, 
the output angles (I,V), and their derivatives. The self-coupling con- 
troller units, #i(s) and B 2 (s), however, depend not only on these plant 
parameters and state variables, but also on the choice of the channel 
input-output transfer functions, Di(s) and D 2 (s). 

Since both B\ (s) and B 2 (s) have the same functional form, we will use 
the subscript i = 1,2 to refer to either channel. From (23), B,(s) is 
written in the simplified form, 

™-(m- l )-i~k' (24) 

where for i = 1: Ci - Ft/go , .A = F i /(g 1 + h), 

and for i = 2: c 2 = F 2 / (g, + ag 4 ), f 2 = F 2 / (g 3 + h 2 ). 

Accuracy, fast response, stability, and a practical feedback structure 
are the general objectives which should be mutually satisfied to the ex- 
tent possible in the choice of the channel transfer function, Di(s). More 
specifically, we consider the following requirements: 

i. Each channel should have no steady-state error for step and ramp 
inputs. 

ii. The feedback synthesis of 5,-(s) should employ signals proportional 
to the output angle, its velocity, and acceleration, but no higher deriva- 
tives. 

iii. The error for a sinusoidal input, having an angular velocity no 
greater than »>, and an acceleration no greater than 7, , should not exceed 
the allowable value, p, . (Appropriate numerical values of v, 7, and p to 
be specified for each channel.) 

iv. The choice of D,(s) should achieve a good compromise between the 
competitive aims of fast transient response and small noise band- 
width. 

To satisfy requirements (i) and (ii), the channel transfer functions 
must be of the form 

dm = , + £Xt + * <25) 

where the coefficient c, is the same as in (24), but &,• and rf, are available 
as design parameters. Using (25) in (24) yields 



^^[(M-^ + £«] 



(26) 
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whore c, and /, are defined below (24) in terms of the linear plant con- 
stants shown in Fig. 10. 

The error transform for either channel is 

^W-ft-IWUW- yJffigk (27) 

where ip r is the channel reference input (either I r or V,). Considering 
requirement (iii), tp r is a sinusoidal signal, 

(p r (t) = <pi(w) sin tot, 
where the peak amplitude is given by 

The steady-state error for this input will not exceed the allowable value, 
Pi , provided 

I 1 - Di(ja) | ^ Pi/Mo)- (29) 

Both sides of (29 ) are shown in the log amplitude-log frequency plot of 
Fig. 11, using (28) and the asymptotic straight-line approximation of 
(27). Since at low frequencies, 

| 1 - D t (ja) | = bj/di, (30) 

it follows that requirement (iii) imposes the design constraint (see Fig. 
ID, 

bi/di^pihi (31) 

which shows that for a given tolerable error, p,- , the limit on the parame- 
ter ratio is imposed by the peak acceleration, 7,- , rather than the peak 
velocity, v { . 

In considering the transient response of each channel, we make use of 
a computer study made by J. F. Kaiser on the step and ramp response 
of the equivalent system 

DM = "r + 4 T , (32) 

s 6 + nw s- + t/cVs + coo 3 

which has a noise bandwidth related to the coefficients by, 

2 

7_ 

!(W - 1) 
From this study, the parameter values which seem to give the best com- 



<B-/ |D(ico)| 2 ff= * + * «■• (33) 

•>-•>» 2 ( vu — 1 ) 
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Fig. 11 — Amplitude-frequency sketch of error functions in relation (29). 



(34) 



promise between good transient response and small noise bandwidth lie 
in the region, 

[2.3 < r? < 2.7 \ 
\1.8 < n < 2.3/ 

within which the noise bandwidth varies from about 0.9 coo to 1.1 oj . 
For small bandwidth, « should be chosen small. However, since bi = nm 
and di = m, the lower bound on u> to satisfy (iii) is, from (31), 

coo 2 ^ m(t,/p,). (35) 

Furthermore, the choice of <o is also linked to the physical system con- 
stants, since c, = t/ojo . 

The detailed choice of the numerical values in the controller design 
cannot be made until all the values of the system constants and the oper- 
ating requirements are known. However, the specific structure of the 
non-interacting controller designs can be given in terms of the channel 
gains F\ and F 2 , the plant parameters go through g 6 (see Fig. 10), and 
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the variable parameters bi , di , 62 and d 2 . This controller design for the 
normal autotrack mode of operation is shown in Fig. 12. The design of 
supplementary control action for special tracking modes is considered in 
the following section. 

4.2 Near Zenith Control Modes 

As noted in Section II, there are two possible tracking modes for a given 
satellite track. For later reference, let "mode 1" (h , Vi) be the tracking 
mode where ^ / ^ 180° and "mode 2" (Z 2 , V 2 ) be the tracking mode 
where 180° ^ I ^ 360°. These two modes are illustrated in Fig. 13 as a 
function of time for the tracking of a satellite in a circular equatorial 
orbit of 0000 miles altitude with the antenna site located at 2° latitude. 
As can be seen from this figure, the maximum vertical axis angular 
velocity required to stay on track, which we shall call y max , occurs at 
the point of maximum azimuth tracking rate and maximum elevation. 
As the maximum elevation angle approaches 90°, F ma x will eventually 
exceed the maximum vertical-axis velocity of the antenna, £. 

Fortunately, there is a factor which reduces vertical axis tracking re- 
quirement. The antenna has an "on-track" beamwidth,* 2£, so that it is 
possible to track without the antenna pointing directly at the target. 
Wo can utilize this "on-track" beamwidth in the following manner. For 
any satellite path which passes within £ degrees of zenith, it is possible 
to switch tracking modes without any interruption in communications. 
This situation is illustrated in Fig. 14, where Vi , V 2 tracking modes 
near one of the switchover points are plotted as a function of time for an 
assumed beamwidth of 0.2 degrees, a maximum satellite elevation angle, 
Em*x = 89.9°, and a satellite altitude of 6000 miles (circular equa- 
torial orbit). The solid lines show the vertical axis angle if the satellite is 
on boresight. The two pairs of dashed lines represent the allowable varia- 
tion in the vertical axis pointing angle due to the beamwidth of the an- 
tenna. If the vertical axis angle is anywhere within the area between 
the dashed lines, communications can be maintained with the satellite. 
Since the tracking areas have a point of intersection at V = V, , a smooth 
transition between tracking modes is possible without a lapse in com- 
munications. For satellite tracks where £ max > 89.9°, the two tracking 
mode areas will have a large area of intersection rather than a single 
point of intersection as in the limiting case discussed above. 

This operation of switching tracking modes will be referred to as the 
"switch-mode." If the switch mode is employed, the vertical axis is to 

* That is, the beamwidth for which the signal-to-noise ratio at the receiver is 
considered sufficient for communication objectives of the system. 
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400 




-0.3 0. 

TIME IN HOURS 



fig. 13 — Antenna pointing angles for satellite in circular equatorial orbit of 
(iOOO mile altitude with antenna station at 2° latitude. 



track keeping the satellite on horesight until, at time t al in Fig. 14, 
the vertical axis angle of the tracking mode equals the angle at which 
switchover from one mode to the other will occur. This switchover angle, 
V, , is equal to the azimuth angle at the maximum elevation point and 
is computed in advance from the predicted satellite orbit. At /„i , a verti- 
cal error signal (V„ — V) is employed to keep the vertical axis angle at 
V M • This error signal (V, — V) is maintained until, at time t,t as in Fig. 
14, the vertical axis angle for horesight tracking of the second tracking 
mode is equal to V, . At this point the switch mode vertical axis error 
(V s — V) is replaced by the horesight tracking error (V r — V) and 
tracking is continued using the second tracking mode. 

Graphical investigation of this switchover process indicates that, de- 
pending on the satellite track, a particular tracking mode should be 
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-0.04 



TIME IN HOURS 



Fig. 14 —Tracking angles of vertical axis near one of the switchover points 
for E mn% = 89.9°, £ = 0.1°, and 6000 mile circular equatorial orbit. 

chosen for tracking from the horizon, where the satellite is first acquired, 
to the maximum elevation point and the other mode from the maximum 
elevation to the horizon. The choice is made by determining, from orbit 
predictions, the initial vertical axis angular velocity for each tracking 
mode when the satellite initially appears at the horizon. If the initial 
vertical axis tracking velocity of mode 1, V H1 , is greater than the initial 
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tracking velocity of mode 2, V in , then mode 1 is used first and the switch- 
over is made to mode 2. The opposite procedure is employed if Vh2 > 
Via ■ Although the tracking velocity requirement is greater at the 
horizon using this procedure, the tracking velocity requirements near 
zenith are reduced. 

For satellite tracks where E max < 90° — £, one cannot switch tracking 
modes without losing communications for a brief period. Therefore, one 
desires to track the satellite by remaining in the same tracking mode. 
However, instead of pointing directly at boresight when tracking the 
satellite, one can again utilize the "on-track beamwidth" to point off 
boresight and still be within the antenna beam as shown in Fig. 15. By 
following the tracking path, V, illustrated in Fig. 15, it is possible to 
significantly reduce the peak vertical axis velocity requirement from the 
velocity required for boresight tracking. This procedure will subsequently 
be referred to as the "slant-through" mode. 

In the design of the antenna drive system, an upper limit is needed on 
the maximum vertical axis angular velocity required to follow the slant- 
through path, V, for all circular satellite orbits of a given altitude, i.e. 
the worst case slant-through speed, 8. Since it is expected that the com- 
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Fig. 15 — "Slant through" path for (iOOO mile altitude, circular polar orbit 
with #„ m x = 89.893°. 
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munication satellites to be tracked will be launched in the same direction 
as the earth's rotation, the worst case will occur when the inclination 
angle, a, equals 90° (a polar circular orbit) and E max = (90° — £). 
(This elevation is chosen since if E max is larger than (90° - £) the switch 
mode will be employed.) An approximate expression for S is derived in 
Appendix C. Shown in Fig. 16 is a graph of £ as a function of satellite 
altitude, x = (r — 1), for antenna half-beamwidths of 0.1° and 0.2°. 
For example, the continuous tracking of satellites in nearly circular 
orbits of 6000 miles altitude with an antenna half-beamwidth of 0.1° 
would require a vertical axis slewing capability of approximately 1 rpm. 
A major factor affecting this capability is the gear ratio used in the 
drive system, the choice of which depends also on other important con- 
siderations such as tracking performance at very low speeds, immunity 
to disturbance torques, minimization of reflected load inertia, and drive 
power requirements. 

This tracking capability should not be difficult to achieve for medium 
and high altitude satellite systems. For low altitude systems, continuous 
tracking may be achieved with the same speed capability by temporarily 
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Fig. 10 — Vertical axis velocity requirement as a function of satellite altitude 
for antenna half-beaniwidths £ = 0.1° and 0.2°. 
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broadening the antenna beam during the near-zenith portion when 
signal strength conditions are favorable. 

The slant-through mode can be performed by biasing the vertical axis 
error signal with another signal (V — V T ) between t x < t < t 2 as shown 
in Fig. 15. V is the computed vertical axis angle which will program the 
"slant-through" path within the antenna beam width. V r is the computed 
vertical axis angles for boresight tracking determined from orbit predic- 
tions. Since V r — V T , the resultant error signal looks like (V — V). 

The above discussion indicates the following tracking strategy. From 
orbit predictions of the satellite path, the elevation values (E) as a func- 
tion of time are scanned. If the maximum elevation, E max , is greater 
than or equal to (90° — £°) the "switch" mode will be employed as ex- 
plained above. If E max < (90° — £), but (V mB . x ) > S, the slant-through 
mode is employed at time t\ in Fig. 15. If neither of these special modes 
is required, normal tracking will be employed. For the normal and slant- 
through modes, it is desirable to use the tracking mode which has the 
smaller vertical axis velocity required when initially acquiring the satel- 
lite at the horizon. The inclined axis velocity requirement at satellite 
acquisition is the same for either mode. If V B i > V m , then mode 1 will 
be used for tracking and the opposite choice will be used if V H \ > V H 2 . 

The mode tracking strategy is shown in flow chart form in Fig. 17. 
Note that no special control signals are needed for the inclined axis. The 
initial inclined and vertical pointing angles when the satellite appears at 
the horizon, I H and V a , can also be determined from orbit predictions 
and (8). 

V. DESIGN RESULTS AND CONCLUSIONS 

The second phase of this design study consists of a simulation and 
design evaluation program on a hybrid analog-digital computer facility.' 
The major objectives of this program are to verify, improve, and if 
possible, simplify the basic controller strategy developed in Section IV. 

The analog computer portion of the facility is being used for the 
simulation of the controller and the antenna dynamics. The experimental 
results discussed here were based on a simulation of an open cassegrain 
antenna with a 56-foot aperture, an overall height of about 70 feet, and 
a total weight of about 100 tons,* using two 25-hp hydraulic motors for 
the vertical axis drive and two 10-hp motors for the inclined axis drive. 

The closed-loop response of the antenna drive system using the basic 
controller design shown in Fig. 12 was tested using step and ramp in- 

* Subsequent design modifications have changed the weight and the antenna 
dynamics somewhat, but not enough to significantly affect the simulation results. 
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Fig. 17 — Flow chart of tracking strategy. 



puts, I T and V r , as well as wind disturbance torque inputs. The design 
for each channel was based on the approximate model transfer function 
(32) with parameter values, 77 = 2.4, p = 1.9, and a range of values of 
o) from 5 to 40 sec. -1 . Because of the low accelerations required for the 
expected satellite tracks, the constraint (35) was easy to satisfy with the 
available range of loop gains; therefore the major considerations in the 
choice of o> were good transient response, steady-state accuracy, and 
immunity to wind disturbances. Satisfactory performance with respect 
to the objectives of zero steady-state error in tracking constant velocity 
inputs and minimum channel interaction was achieved in the experimen- 
tal design with state feedback gain adjustments close to the nominal 
values computed for the expressions in Fig. 12. 

The experimental data obtained from this design study has established 
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the feasability of operating the open cassegrain antenna under severe 
wind conditions without a radome. A series of tests using simulated 
wind loads corresponding to a 40-mph gale with gusts exceeding 80 mph 
have indicated that the control system is capable of maintaining the 
antenna beam on-track with both a mean and rms error less than 0.002 
degrees, which is about l/100th of the nominal beamwidth of the an- 
tenna. 

The future test program in this design study will employ the digital 
computer portion of the hybrid computer facility to simulate satellite 
tracking data, autotrack error detector signals, and the necessary co- 
ordinate conversions for resolving the error signal inputs and the angu- 
lar outputs of the non -orthogonal axes of motion. This will provide a 
complete simulation of the overall autotrack system shown in Fig. 2, 
and will allow testing and evaluation of the overall control strategy 
for non-orthogonal mounts, including the near zenith control modes 
discussed in Section 4.2. 
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APPENDIX A 

We consider the polar coordinates [r(t),d(t)] of the satellite in the 
orbital plane to be given, and then make the sequence of transformations 
necessary to relate the tracking coordinates (A,E) to these orbit parame- 
ters: 

A.l Transformation from orbit plane to equatorial plane 

We define rectangular coordinates XYZ and X'Y'Z' as shown in Fig. 
18, with the X and X' axes coincident with the line formed by the inter- 
section of the orbit and equatorial planes. The satellite moves into the 
northern hemisphere at the + X-axis. We define the following (see Fig. 
18): 

r(t) = range from center of earth to satellite at time t. 
8(t) = angle which r(t) makes at time t, measured from the X'-axis. 
a = inclination angle of orbit plane with equatorial plane.* 



* Should not be confused with antenna incline-angle, a, denned in Section II, 
since they will not be used in the same context. 
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Z 

z 




Fig. 18 — Satellite coordinate conversion — Orbit plane to equatorial plane. 



\J/ S (l) = longitude angle of satellite eastward from X-axis. 
\ t (t) = latitude angle of satellite northward from equator. 

The transformation from the orbit plane to equatorial plane corresponds 
to a rotation about the X-axis through an angle a, so that 



M 


r 


Y 


= 


Z 





1 

cos a —sin a 
sin a cos a 



(X') 
Y' 
Z' 



(36) 



Expressing the rectangular coordinates in the equivalent spherical polar 
coordinates, we have 



r cos X, cos yf/ s 

r cos \, sin yf/ s 

r sin A, 



1 

cos a —sin a 
sin a cos a 





r cos 6 | 




r sin 6 




l 



(37) 
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A.2 Transformation from equatorial plane to horizon plane at Antenna 
Station 

The rectangular coordinates (xyz) are located with the origin at the 
antenna site as shown in Fig. 19, with the +.t--axis pointing northward. 
The notation used here is as follows : 

R = radius of the earth (assumed constant ) 

i> = longitude angle of antenna site, measured eastward from X-axis 

X = latitude angle of antenna site, measured northward from equa- 
tor 

p = slant range to satellite from anetnna site 
A = azimuth angle measured CW from .r-axis (North ) 
E = elevation angle measured up from horizon (.r?/-plane ) 

The transformation from XYZ coordinates to xyz coordinates can be 




EQUATOR 



Fig. 19 — Diagram for the transformation of coordinates from the earth- 
centered system (XYZ) to the local antenna site coordinates {xyz) and pointing 
angles (A,E). 
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written as a sequence of two rotations and a translation, as shown in 
Fig. 19. The first is a CW rotation about the Z-axis through an angle 
(180° — \J/), which can be written as 

f— cos \J/ —sin \f/ 0] 



sin ^ — cos ^ 











(38) 



The second is a C W rotation about the F-axis through an angle (90° — X ) 
as indicated in Fig. 19 : 



sin X cos X 
1 

— cos X sin X 



r 



Y 



2 



(39) 



Thirdly, a simple translation along +z'-axis a distance R gives, 



X 




X 




'°1 


y 


= 


1 

V 


- 





z\ 




, z J 




R. 





X 




0] 


= 


1 

y 


- 







i 

< z 




R. 



or, in terms of the spherical coordinates at the antenna site 
p cos E cos A 

— p cos E sin A 

p sin E 
Finally, combining (37) through (40), we have 

p cos E cos A = — r (sin X cosxf/ cos 6 + cos a sin X sin $ sin0 
— sin a cos X sin 6) 
— pcosE sin A = r (sin \f/ cos 6 — cos a cos ^ sin 0) 

p sin E = r (cos X cos ^ cos 8 + cos a cos X sin \f/ sin 
+ sin a sin X sin 0) — R. 
Separating these variables gives the desired transformations : 



(40) 



(41) 



AUTOTKACK CONTROL SYSTEMS 



1397 



Slant range to satellite: 
P = [r 2 + R 2 — 2rR (cos \ cos \p cos 8 + cos a cos X sin \f/ sin 5 

+ sin a sin X sin 0)]* 

Azimuth angle to satellite: 

. _i / sin ^ cos — cos a cos \1/ sin \ 

A = tan ' ^ 



sin X cos iA cos 8 + cos a sin X sin ^ sin 
— sin a cos X sin 8, 



(42) 



(43) 



Elevation Angle to satellite: 



E = sin ' < 



r(cos X cos ^ cos + cos a cos X sin \f/ sin 

+ sin a sin X sin 8) — R 



[r 2 + /i! 2 - 2rfl(eos X cos iA cos 

-f- cos a cos X sin \f/ sin + sin a sin X sin 0)]' 



(44) 



where -tt/2 g E ^ tt/2. 

In these expressions the constant factors are the orbit inclination angle, 
a, the antenna site latitude, X, and the earth radius R. The antenna site 
longitude angle, \f/, can be considered constant if the effect of the earth's 
rotation is neglected, but more generally it will have the form 



iKO = tAo + fi/ 



(45) 



where J2 = earth's angular velocity = (ir/12) rad/hr. To simplify the ex- 
pressions we shall measure distance in units of earth radius (e.r.) so that 

(46) 



R = 1 e.r. 

Further, to emphasize the constants, let 

C a = cos a, C\ = cos x\ 
S a = sin a, S\ = sin Xj 



(47) 



Using the notation of (45), (4(i) and (47), the tracking angles can be 
written 



A = tan 



/ sin \p cos 8 — C a cos \f/ sin 8 1 

[«S'x(cos yp cos 8 -\- C a sin \p sin 8) — S a C\ sin df 



E = sin" 



|?'(Cx(cos \f/ cos 8 -f- C a sin \f/ sin 8) 

+ & a S x sin 0) - 1 



[I + r 2 - 2r(C x (cos tA cos 

+ ('„ sin \P sin 0) + S a( S x sin $)]! 



(48) 



(49) 
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where r, 9, and ^ arc in general varying with time. The denominator in 
(49) is the range from antenna to satellite in units of earth radii. The 
numerator in (49) gives the necessary condition for the satellite to be 
above the horizon at the antenna site, namely : 



(50) 



r[Cx (cos ^ cos + C a sin ^ sin 0) + S a Sx sin 0] ^ 1 => E > 

where, of course, r is the radius vector of the satellite in units of earth 
radii. 

During the period when (50) is satisfied, the A (t) and E(t) given 
above, as well as their time derivatives, give the required information 
to evaluate the antenna drive angles and rates, from (8) and (10) in Sec- 
tion II. For non-circular orbits, the expressions for (r,0) as functions of 
time must still be determined, and in general a computer routine 4 would 
be used. The main usefulness of the analytical expressions (48) and 
(49) is in estimating tracking requirements for particular orbits where 
complete data are not needed. 

A.3 Tracking angles and rates for circular orbits 

For a given inclination angle, a, of the circular orbit plane (see Fig. 18 ), 
and a given latitude, X, of the antenna site, the tracking angles A (t) and 
E(t) can be obtained from (48) and (49) as explicit functions of time by 
substituting the time functions </> (t) , given in (45), and 

B = ut + O (51 ) 

where ( , = 0(0 at arbitrary time reference, t = 0, and the constant 
angular velocity of the satellite is given by 

u = kr~ m (52) 



where 



or 



r « l -l h/R, 

h = satellite altitude 

R = earth radius 

A- = (g/R)' 3 

g = accel. due to gravity 



k = 4.47 hr _1 = 1.24 X 10 3 sec '. 



The particular case of the circular equatorial orbit yields the simplest 
expressions for azimuth and elevation angles and rates. Letting a = 0, 
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we obtain from (48) and (49) the angle expressions: 



i tan <p (t) 



[A (0]_o - A*(t) = tan" 1 - ^^ (53) 

{ sin X J 

rn/Ai n*/\n • -i ( r COS X COS a> (<) — 1 \ /_-n 

[E(t)] a=0 = E*(t) = sin 1— —r^ (54) 

([1 4- r — 2r cos X cos tp{t)] *J 

where 

<p(0 = (co - 12)£ + 0o - i/'o . (55) 

Taking the time derivative of (53), the azimuth rate for circular equa- 
torial orbits is 

[A(t)U m A*(t) = . 2 ... ~? n * r-r^ (56) 

sin 2 <p(t) + sm 2 X cos 2 <p{t) 

where, from (55), <p = co — U = constant relative angular velocity of 
the satellite with respect to the earth. The maximum azimuth rate 
occurs when <p(t) = 0, and is given byf 

I 4„„* I = ^» (57) 

sm X 

where, from (45) and (52), 

<p= ( w - Q) = 4.47 r~ m - (tt/12) rad/hr (58) 

the maximum elevation, E, also occurs when <p(t) = 0, and has the 

value, 

rcosX-1 (ri(|) 



[1 4- r 2 - 2r cos X]*j ' 

Differentiating (54), the elevation rate for the case of circular equa- 
torial orbits is, 

jpm( f \ _ — r cos X sin <p(t) (r — cos X cos <p(t))<p , . 

W (1 - cos 2 X cos 2 ^(0) Hi + r 2 - 2rcosXcos^(0) ' 

where <p(t) and and <p are given by (55) and (58), respectively. We 
are interested in the maximum value of E when the satellite is visible, 
i.e. when r cos X cos (pit) > 1. Once the parameters r and X have been 
specified, this can be determined from (60). 



t If / = 180°, then V max = A max . 
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APPENDIX B 

We derive the inclined and vertical channel error signals as functions 
of the waveguide error signals and the controlled antenna angles (I,V) 
in this appendix. The desired error signals are 

6! A I r - I 

(61) 
e 2 = V r - V . 

Since we desire to make ei(t) to be within some small specified toler- 
ance, we assume 

Bini(/ r -/) = «i/2 

and (62) 

sin \(I T — I) — sin /. 

Using (62) and the trigonometric identity, 

sin a (A + B) sin \{A - B) = -|(cos A - cos B), 

we can write 

(cos / — cos I T ) 



*i = 



sin / 

Using (8) and (63), one obtains 

(sin E r — sin E) 



«1 = 



(63) 



(64) 



b 2 sin J 

Using (1), the trigonometric identity for the difference between sin E r 
and sin E, and assumptions similar to those used above 

_ (E r — E) cos E _ e v cos E ,__v 

61 6 2 sin / b 2 K v sin / ' 

From (7) and some algebraic manipulation, we can write (65) as 

(66) 



ei = 



cD + '""©]■ 



The vertical-axis error signal, e 2 is found using a similar procedure. 
From (61) and (8), 

€2 = (V r -V) = A r - A(/ r ) - A + A(J). (67) 

From (1), we have 



AUTOTKAG'K CONTROL SYSTEMS 1401 



e 2 = „ 6h „ + A(7) - A(7 r ). (68) 



K h cos E 
Using the trigonometric identity, 



sin (A (Ir) -MD) 



tan A(/ r ) - tan A (7) = 7,7 =^ff (69) 

cos A(7 r ) cos A(7) 

and, assuming that for normal tracking (A(/ r ) — A (7)) is small, we 
obtain 

A(/ r ) - A (7) = |[tan A(7 r ) - tan A (7)] 

X Jcos[A(7 r ) + A(7)] (70) 

+ cos[A(7 r ) - A (7)]} 

using (69) and the trigonometric identity for the product of two cosines. 
Using (7) and assuming that 



~ \ [ tan © = tan (j). 



(71) 



a[\ + cos 1] 
cos[A(7 r ) - A (7)] = 1 
cos [A(7 r ) + A (7)] = cos2A(7) 
(70) becomes 

A(7 r ) - A(7) = - | [1 + cos 2A(7)][1 + cos if. (72) 

After some additional algebraic manipulation using (7) and trigonomet- 
ric identities, one obtains 



A(7,) - A(7) = -ae,[l - sin 7?] \ (73) 

Using (68) and (7), (73) becomes 



r 



e, =* S- 4 b 2 sin 2 7(l + a 2 cot 2 £ ) I + ^ esc 2 ( £ ) . (74) 



APPENDIX C 



An approximate expression for S is derived in this appendix. Refer- 
ring to Fig. 20, one can say 

sin [(,r/2) - £] = [r cos <fo - l][(r cos ft, - I) 2 + (r sin ^ )T* (75) 
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BEAMWIDTH 



SATELLITE 



Fig. 20 — Tracking geometry for worst case slewing speed calculation. 
Since £ and \f/ are small angles, (75) can be written as 



-f- 



1 - 



r\pQ 



1 + 



r\po 



2(r - 1)_ (r - 1) 2 _ 



Since r\p 2 / (r — 1 ) 2 <5C 1, then (76) can be expressed in the form 



r = 



rypi 



(r - 1. 



+ 



ripa 



+ 



(76) 



(77) 



where the remaining terms of the series are negligible compared to the 
first two terms in (77). Solving (77) for £, 



or 



t = 



lAo — 



r^t 



(r - 1) 
(r - 1) 



t 



(78) 



Using (48), (49), (8) and (78) for the special case of a circular 
polar orbit of 6000 mile altitude (r = 2.5 e.r.) and £° = 0.1°, the vertical 
axis angle pointing directly at the satellite and boundaries on either side 
of the vertical axis angle which are allowable for tracking of the satellite 
were graphically determined and plotted in Fig. 15. The maximum ve- 
locity requirement, o, was graphically determined by picking the 
minimum slope as shown in Fig. 15. For the particular example plotted, 



S = 1 R.P.M. 



(79) 
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The velocity requirement for "boresight tracking," S, is obtained by 
taking the time derivative of (48) with S a = C\ = 1 and evaluating 
at / = 0. The result is 

S = | y mas u_9o° = «/*<> (80) 

for small ^ angles. Using equation (78), one can write 

Using (81) for a 0000 mile polar circular orbit, the comparable slewing 
requirement for "boresight tracking" is 2.85 rpm. Graphical plots for 
other orbit inclinations, altitudes, and antenna beamwidth angles indi- 
cate that tracking on the "edge of the antenna beamwidth," as illus- 
trated in Fig. 15, reduces the required slewing capability for a given 
circular orbit by approximately g. Therefore, an approximate expression 
for § is 

S = jg. (82) 

Using (81) and (52), (82) can be expressed as 

S . ™Mg, (83) 

r*(r — 1)£ 

where r is in earth radii and £°, the half-beamwidth angle, is in degrees. 
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